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ABSTRACT 

During deep-drawing, two issues manifest themselves 
that are due to normal and planar anisotropy in the sheet. These 
are a dependence of the Limiting Draw Ratio (LDR) on the 
average thinning ratio (i?-value), and a dependence of ear 
formation and thinning around the circumference on the 
variation of R-value in the plane of the sheet. The quadratic 
(1948 Hill) yield surface has been applied to these issues and it 
has been demonstrated that there are numerous higher exponent 
yield criteria that may more closely duplicate experimental 
trends. These predict varying degrees of R- value dependencies 1 
of uniaxial yields and strength ratios in multiaxial loading paths. 
The result of this is that there are sometimes subtle and 
sometimes substantial differences in the predictions of the 
various yield surfaces on deep drawing regarding LDR and 
earing. Additional differences arise due to the way the shear 
term (45-degree yield) and the 0-degree vs. 90-degree strengths 
are treated in each criterion. These dependencies (in-plane 
strengths generated by the yield surfaces) are shown to affect 
the results of both LDR and earing during cupping. In 
particular, the 1979 Hosford and 1989 Barlat (Tricomponent) 
criteria, although identical for die normal anisotropy (planar 
isotropy) case, are strikingly different for cases where there are 
differences at 45 degrees and 90 degrees from the rolling 
direction. 



INTRODUCTION AND YIELD SURFACES 

In recent previous works (Logan, 1995 and 1996) we 
have outlined the implementation and use of various anisotropic 
yield criteria in explicit dynamics finite-element codes such as 
DYNA3D (Whirley, 1991) or ALE3D (Couch, 1995). Those 
yield criteria have included quadratic (Hill, 1948), higher- 
exponent (Hosford, 1979), and TriComponent (Barlat, 1989). 
All three of these yield criteria can be expressed in terms of a 
uniaxial flow stress X together with Strain ratios from uniaxial 
tension or compression tests in various directions. We postulate 
the material coordinate system where a=rolling direction, 
&=transverse direction, and c=normal direction to the sheet, 
with strain ratios: 



R=R 0 = e b /e c 

= R 45 = E „/ £ c 

P=R 90 = eJe c ...(1) 

Qca = 

Here, the a'-b' coordinate system is obtained by rotating the a-b 
system 45° counterclockwise, the b"-c” system by rotating the 
b-c system 45° counterclockwise, and the system by 

rotating the a-c system 45° counterclockwise. R is measured in 
an a-direction uniaxial test, Qab in an a -direction test, P in a b~ 
direction test, Qbc in a (hypothetical) fr' '-direction test, and Q ca 
in a (hypothetical) a" -direction test. In reality, Qbc and Qca are 
not overly critical to sheet forming except with very tight 
bending radii, and a von Mises approximation is often used as it 
is here. Having measured or approximated the above strain 
ratios, we have all the needed constants for the quadratic 1948 
Hill (H48) yield criterion: 

' —2 F(cr b -<T c f+Gi<T c -<7 a ) 2 + H (<T a -<T b ) 2 +D 

a - — ...( 2 ) 

Eqn. (2) relates the effective stress to the three normal 
components of Cauchy stress, with the term D containing die 
shear stress terms: 

D^lLai+lMai+lNai ...(3) 

The values for the constants in Eqns. (2) and (3) can be 
expressed in terms of above in-plane strain ratios for sheet 
forming, with the following additional relations needed: 

F-S-R/P 

G-l 

H=R 
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L = (Q bc +i)(R + l) 

M=(Q ca +i)(R + S) 

N = (Q ab +i)(l + S) 

-(5) 

The quadratic 1948 Hill criterion permits the relatively 
simple calculation of the ratio, X{6)/X, of the yield stress in a 
direction at an angle in the plane of the sheet to the rolling 
direction, as well as the calculated F-value in that direction: 

X(6) 

X 
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many of our implementations refer to 'm' as a strain-rate 
exponent Note also that to avoid confusion we have expressed 
the coefficient (2-c) explicitly: 

(2 - c){\K, + K-X +|^,-A' 2 f}+ c\2KX 
G ~ 2 
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c = cos(0) 



s — sin(0) 

The 2?-value relation giving R(Q) is used in the 
implementation of the 1979 Hosford equation, extending 1948 
Hill to a non-quadratic form with values of the exponent a in 
the range of a= 8 for fee (Hosford, 1979), and a =6 for bcc 
(Logan, 1980) metals: 



The Cauchy stresses must be defined to allow for a third 
(normal) stress, even though the implementation here is for the 
plane-stress shell element. This is because the plane-stress 
material routine is iterative so that even though the normal stress 
vanishes at convergence, we must recognize its presence during 
the iterations. Parameters c, h, and p* may be defined in the 
current notation as follows: 



_ F(<j 2 - (j 3 ) a + G{a 3 - <7,)° + H{a { - a 2 ) a 
R + l 



( 8 ) 
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Eqn. (8) is the second yield surface to be used in this 
study, with its implementation in DYNA3D discussed earlier 
(Logan, 1995). The advantages of the higher exponent in 
eqn. (8) have been discussed at length, and some of them are 
reviewed below; the higher exponent in itself allows an 
adequate fit to many (though not all) measured yield surfaces. 
However, difficulties arise in the implementation of the 1979 
Hosford criterion for cases in other than principal stress/strain 
space. This stems from the lack of shear terms in the criterion. 
Eqn. (8) must remain in principal stress space to be used 
without spurious results and non-convexity problems. To do 
so, we must make a key approximation in updating the Cauchy 
stress tensor when using _eqn. (8). The key necessity is a 
rotation to the principal stress coordinate space. It is 
straightforward to find and rotate the stress tensor to principal 
stress space, but to express eqn. (8) in that space we need to 
rotate the constants F, G, and H as well. To do so, we assume 

the dependence R(Q) as given for 1948 Hill by eqn. (7). The 
principal stress direction does not noxmally coincide with either 
the material (rolling and transverse direction) coordinate 
system, nor with the axes of principal strain. We assume that 
the axes of principal stress and strain approximately coincide, 
although for planar isotropy we know that they normally will 
not. However, this assumption, which leads us to ignore cross- 
terms in the constitutive matrix, is believed to lead only to small 
errors for the degree of anisotropy observed in most sheet 
metals. This is the assumption we use to update the stresses for 
the 1979 Hosford criterion. 

In order to circumvent the stress space limitation of the 
1979 Hosford criterion, Barlat and Lian (Barlat, 1989) 
introduced a criterion which offers all the advantages of the 

1979 Hosford for the case of normal anisotropy (AF=0) but 
permits the introduction of a coupled shear term while retaining 
convexity of the yield surface and coordinate system invariance. 
This criterion is expressed below as in (Barlat, 1989), except 
that we retain the use of 'a' as the yield criteria exponent since 




R(l + P) 

P(l + F) 
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V (1 + XX2+C) 



...( 10 ) 



The value of p is needed for the shear term in eqn (10). 
In the case of a- 2, we have p=p*. However, this is also the 
case where the criterion reduces identically to 1948 Hill and is 
thus of interest only for verification. In general, the value of p* 
must be found iteratively as described in (Barlat, 1989). 
However, this has not been found to be a drawback for the case 
considered thus far where these coefficients are constant 

Given either of the choices of yield surface, the next 
step in implementation involves the stress update. For the 
isotropic case, this is conveniently done using the radial return 
method as discussed by Krieg and Key (Krieg, 1976). 
However, this method cannot be used directly for anisotropic 
plasticity. Thus, the incremental method described in (Bathe, 
1982) and others is used in the DYNA3D and ALE3D 
implementations. Due to the small strain increments typical in 
an explicit dynamics code, this again is not a major drawback in 
computation speed. 

If we consider the three yield criteria described above, 
that is the Hill 1948 (H48), Hosford 1979 (H79), and 
Barlat/Lian 1989 (B89), we can gain a quick insight into the 
tendencies of each by examining some of the simpler yet 
important stress ratios calculated by each. For example, at 45 
degrees to the rolling direction, H48 predicts: 
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Similarly, H79 predicts a similar and simple relation at 
45 degrees to the rolling direction : 



W _ X(45°) 
X~ X 



R + l 

CL+i 



1 % 



(VR = P) 



...(12a) 



Butler (Wilson, 1962); the contrived data sets have been 
adjusted to show particular numerical points. The numbers 
chosen are representative of available sheet materials but 
adjusted to provide exactly R=P in one case (material Sb), and 

exactly AR=0 in another case (material Nm). This allows us to 
clearly separate the subtleties of the different yield criteria. 



For the B89 criterion, there is no explicit formula for the 
ratio ( W/X) ; its value tends to be intermediate between that 
predicted by H48 and H79. The value may be approximated by: 



W = X(45°) 
X X 



R + l 

0* + i 




(VjR = P) 



...(12b) 



We will see the manifestation of these trends below in 
the earing results. Since the exponent ’a’ is still a= 2, Eqn. (11) 
shows a high dependence of yield stress on orientation in the 
plane of the sheet. This dependence is a likely factor in 
overprediction of earing as in (Logan, 1995), as the strong 
material in the 45 degree direction tends to pull in to form the 
wall of the punch (forming a deep trough), while compressing 
the 0 degree and 90 degree walls (forming high ears). The H79 
equation (12) predicts much milder dependencies of strength 
ratios in stress states and directions other than tension in the 
rolling direction. At 90 degrees to the rolling direction, a similar 
situation exists for general values of unequal R and P: 



Y = X(90°) 
X~ X 



P(R+ 1) 
R(P + 1) 



1 % 



...(13) 



Eqn. (13) interestingly applies to all three criterion, but 
with a subtle twist. It predicts the value of ( Y/X) for H48 with 
a= 2, with a fairly strong dependence. The value of ( Y/X) is 
also predicted for the H79 criterion with o> 2 (typically a=6 or 
a=8), giving a mild dependence of (Y/X) on (R/P). For the B89 
criterion, we can still use eqn. (13), but we must use it with 
o=2 regardless of the exponent chosen for the yield criterion 
itself. In other words, B89 gives the same strong (Y/X) 
variation as H48. This leads to subtle differences in values of 
Limiting Draw Ratio (LDR) and Goss-texture (Malin, 1993) 
earing magnitudes (ears at 0 degrees relative to ears at 90 
degrees). 

In biaxial tension, where o a = a b = B , both 1948 Hill 
and 1979 Hosford predict the relatively simple value: 



B_ 

X 




...(14) 



Clearly, values of a > 2 will give a much milder 
dependence often observed experimentally and thus should 
provide better correlation with LDR as shown in previous 
works. The value of (B/X) for the B89 criterion has a more 
complex form: 



X 



(2~c)(l + h a ) + c(l~~h) a 



...(15) 



Eqn. (15) simplifies exactly to eqn. (14) for the case of 
R-P , thus providing B89 with the same mild dependence of 
(B/X) on R as shown by H79. Graphical examples of these 
trends in stress ratio are shown in Fig. 1 for two materials with 
hypothetical sets of R - values taken from Table 1 below. The 
experimental data sets in Table 1 are taken from Wilson and 
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FIG. 1. COMPARISON OF STRENGTH RATIOS 
PREDICTED USING PLANE STRESS YIELD CRITERIA 
(H48, H79, B89) FOR MATERIALS WHERE R=P (SB), 

AND WHERE R<P WITH A/?=0 (NM). 



Note that in Fig. 1, the H48 criterion predicts strong 
strength ratio dependencies in all cases, and the H79 criterion 
predicts very mild dependencies. In general, the B89 criterion 
agrees with H79 for the B/X ratio, with H48 for the Y/X ratio, 
and is intermediate for the W/X ratio. Note that for the “Nm 11 
material, B89 is again somewhat intermediate between H48 and 
H79 since R«P for this case ; note again the values of R for 
the two contrived materials in Fig. 1, the ’’Sb” and "Nm” sets. 

In the following sections we will demonstrate the effect 
of the chosen yield surface (H48, H79, or B89) on die extent of 
earing in cupping with a "cross-rolled” texture, the nature of the 
LDR dependence for that case and cases where R P, and on 
the earing trends that might be observed for materials fitted to 
these yield surfaces in cases where R±P, 




TABLE L PLASTIC ANISOTROPY FOR CONTRIVED AND 
REAL MATERIAL DATA SETS USED IN THIS WORK . 

Matl \ar 1a P AP j£ 

P Qdb P R R R P P 




EARING: CROSS-ROLLED ASSUMPTION 

To compare the different yield criteria in DYNA3D and 
ALE3D with earing observed experimentally, we compare the 
data obtained by Wilson and Butler (Wilson, 1962) with that 
obtained using simulations with either a =2 in eqn. (2) (H48), 
a =8 in eqn, (8) (1979 Hosford), or a =8 in eqn. (9) (1989 
Barlat). Numerous runs with R = P were made as outlined in 
detail previously (Logan, 1996) for the three criteria outlined 
here, using a 100mm punch and 200mm diameter blank, and 
the R- values as stated in Table 1 under data for Fig. 2. In all 
cases, a rolling-direction yield of 200MPa, tangent modulus of 
250MPa, sheet thickness of t=0.8Gmm, and friction coefficient 
of 0.100 was used with punch nose and die lip radii of 
r =6. 3 5 mm, providing an (r/t)=8 in order to minimize the 
complicating effect of bending. 

The stress ratio dependencies illustrated in Fig. 1 can 
help to intuitively explain the correlation of the various yield 
surfaces with experimental (W&B) earing results shown in 
Fig. 2. As explained above in reference to eqns. (11-12), the 
high (relative to most experimental observations) strength ratios 
predicted by H48 lead to an overprediction of earing on planar 
anisotropy. In this context, we quantify the 4-fold earing with 
ears at 0 and 90 degrees and troughs at 45 degrees as: 

£45 = Ah/h 

AIt/h _ 2*(h0Q-2*h45 + h90) 

(hOO + 2 * M5 + h90) ...(16) 

AR = (£ — 2 * go* + P) / 2 

£=(£ + 2*0, + P)/4 

The quantities £45 and AR/ R are then correlated to 
produce the earing trend lines shown in Fig. 2. As expected 
from its intermediate stress ratio behavior, the earing trend line 
for B89 is intermediate between H48 and H79, the latter of 
which correlates best to the experimental data. 




I ^ OOOO OOOO *-« 

II till 

AR/R 

FIG. 2. PLOT OF EARING (AH/H VS. AR/R) INCLUDING 
WILSON AND BUTLER DATA COMPARED TO ALE3D 
AND DYNA3D RESULTS WITH A =2 (48 HILL), A =8 (79 
HOSFORD) AND A =8 (89 BARLAT) . 



(PLANAR) ANISOTROPY AND THE LDR 

Another strong effect of anisotropy in deep drawing 
occurs with the observed dependence of the Limiting Draw 
Ratio (LDR) on the R-value(s) of the sheet The LDR is defined 
as the ratio of maximum (limiting) blank diameter to punch 
diameter in cylindrical cupping. This dependence was shown 
with experimental data and axisymmetric finite-element 
calculations in previous work (Logan, 1987). In Fig. 3, the 
experimental data in (Logan, 1987) is reproduced, along with 
calculations performed with DYNA3D using the H48, H79, 
and B89 criteria. For this case, with R-Q ab - P— R set 
intentionally, the H79 and B89 should in fact give identical 
results. In fact, there is a slight difference between the two 
which we attribute to perhaps the numerical rounding associated 
with the coincident principal stress and principal strain 
assumption. Nonetheless, both show a much milder 
dependence than H48. Although these trends are essentially the 
same as those shown in (Logan, 1987), the three-dimensional 
implementation now allows us to examine the LDR in a deeper 
way. The experiments performed by Meuleman and reported in 
(Logan, 1987) were divided into three categories, depending on 
the magnitude of the quantity A R/R, as shown in Fig. 3. 
Somewhat surprisingly, there is no discernible affect of AR/ R 
on the LDR as viewed in Fig. 3. To see if the numerical 
simulations and chosen yield criteria concur, we repeated the 
LDR simulations and plotted data points at R = 0.5,/? = 2.0 as 
above with an imposed value of AR/ R = 0.4 to explore the 
affect of planar anisotropy as predicted by each yield surface. 
The H48 criterion shows a significant increase in LDR with 
A R/R = 0.4, in contrast to the experimental data and perhaps 
our intuition as well. The H79 criterion shows only a very 
slight increase in LDR, while B89 actually shows a slight 
decrease. The divergence of H79 and B89 as A R/R is 
increased is not surprising due to their different nature as 
explained in eqns. (12-13) above. What is perhaps most 















striking is the minimal affect of AR/ R on the LDR in general, 
as observed both experimentally and by the two high-exponent 
yield criteria. We must bear in mind that nonzero AR/R will 
inevitably lead to ears and troughs thus lowering the trimmed 
cup height as explored in (Logan, 1996). 




FIG. 3. EFFECT OF R AND AR/R ON THE LDR AS 
SHOWN EXPERIMENTALLY AND BY THE THREE 
YIELD CRITERIA (H48 LABELED A=2, H79 LABELED 
A=8, AND B89 LABELED A =9; EXPONENT IS STILL 
A=8). 



EARING: EFFECT OF UNEQUAL R AND P 

Although the numerical trend lines of earing (E45 in 
eqn. 16) shown in Fig. 2 were obtained with /?=I\ this is rarely 
exactly the case. In order to explore the effect of unequal R and 
P, some of the actual values reported in (Wilson, 1962) along 
with some contrived cases with R<P were run with DYNA3D, 
and the results plotted in Fig. 4a-b. With the simulation results 
connected point-by-point in Fig. 4a, we see that although the 

value of AP = R - P can locally affect the value of Ah/h (E45), 
overall the regression lines shown in Fig. 4b agree with the 
observations of Fig. 2 where R=P . In die cases of unequal R 
and P , we now have another type of earing as well, that in its 
extreme will lead to Goss-texture behavior with ears at 0 and 
180 degrees and troughs at 90 and 270 degrees (Malin, 1993). 
This earing may be quantified by: 



E90 — AA(Q — 90) / h 



A/z(0 - 90) / h 



4(M)Q-/i90) 
(M)0 + 2*M5 + ft90) 



AP=R-P 

P = 0R + P)/2 



-..(17) 



In Fig. 5a-b, we attempt to correlate E90 with the 
quantity AP/ R. Two things are readily evident. First, the 
correlation in general (as noted by the regression values of r 2 ) 
are not very good. Second, the highest dependence is shown by 
B89, rather than the usual winner, H48, in such dependencies. 
Overall H79 shows the least dependence of E90 on AP/ R; 



this at first seems intuitively satisfying, but a simple "flange in 
compression" earing analysis will show the opposite trend, as 
discussed below. In Fig. 6a-b, we correlate E90 instead with 
the quantity AP/ P as defined by eqn. (17) above. There is in 
general still considerable scatter, although even the point-to- 
point plot in Fig. 6a looks smooth for H79. The regression 
lines in Fig. 6b show the same qualitative results as in Fig. 5b, 
with the r 2 coefficients much improved for B89 and H79. 
Alternately, we can correlate E90 with the ratio R/P as shown in 
Fig. 7a-b. Again, the point-to-point plot in Fig. 7a is smooth 
only for the H79 case, and the regression fits are quite good for 
both B89 and H79. The trend of strength of dependencies, B89 
strongest and H79 weakest, remains as it did in Fig. 6. 



DISCUSSION OF CALCULATED E90 TRENDS 

Intuitively, we might at first be satisfied with the mild 
dependence shown by H79 in Figs. 4-7, and surprised by the 
strong dependence of E90 on AP/ P etc. shown by B89. 
Unfortunately, more experimental data is desperately needed in 
this area and will be a topic of future work. Meanwhile, in an 
attempt to numerically resolve these observations, we consider 
here a very simple earing analysis. We will choose two of the 
contrived materials (Nf and Nj) from Table 1, and examine the 
DYNA3D predictions of E45 and E90 in Table 2: 



TABLE 2. TRENDS IN E45 AND E90 EARS PREDICTED 
BY DYNA3D. 



Mad. 


R 


Qab 


P 


Earing 


H48" 


H79 


B89 


Nf 


4 


1 


4 


E45= 


.396 


.259 


.300 


Nj 


1 


2 


3 


E45= 


-.110 


-.060 


-.070 


Nj_ 


1 


2 


3 


E90= 


.091 


.046 


.128 



Material Nf, with AR / R = 1.20, shows die expected trend of 
E45 shown in Fig. 2. That is, H48 shows the most E45, B89 
slightly less, and H79 least of all and closest to the experimental 
trend. The E90 earing is null in each case since R=P exactly. 
The trends for material Nj are particularly interesting, and are 
plotted graphically in Fig. 8 with the value AP/ P = 1.0 for the 
Nj material; the material Nm shows the same trend with its 
AP/ P = 0.5 although to a lesser extent than for material Nj. 
These materials have AR / R = 0.00, hut they do exhibit some 
degree of 4-fold E45 earing, with ears at about 45 degrees and 
"troughs" at 0 and 90 degrees, although in most cases for Nm 
and Nj, the height at 0 degrees is nearly as high as the height at 
45 degrees. Nevertheless, when defined as in eqn. (16), E45 
again follows the familiar trend of H48 highest, with B89 
intermediate and H79 the smallest magnitude. 

Since R«P for material Nj, we can examine the trends 
of E90 with the various yield criteria as shown above in 
Figs. 4-7. For both materials Nm and Nj, the same trends in 
E90 predictions of the three yield criteria are observed; these 
trends continue in the other materials where R^P as well. The 
tendency is that H48 shows an E90>0, nearly proportional to 
A P/P. The B89 criterion shows an even greater E90, instead 
of the reduction vs. H48 observed in predicting E45. The H79 
criterion shows much less E90 than either. From our experience 
with the trend of these three yield criteria in predicting E45, 
B89 would seem to be the outlier if the same ordering (H48 
most earing, H79 least) held true for E90. Our intuition for E45 
was based simply on the predicted value of W/X for each yield 
criterion. However, for E90, the situation is complicated since 
Y ^ X . Thus, our intuition must now encompass the strength 
ratios in Table 3 for material Nj. 







FIG. 7A) E90 CORRELATION, RAW DATA SETS. 



0.5 





FIG. 4B) E45 CORRELATION, REGRESSION FITS. 




FIG. 5B) E90 CORRELATION, REGRESSION FITS. 




FIG. 6B) E90 CORRELATION, REGRESSION FITS. 




FIG. 7B) E90 CORRELATION, REGRESSION FITS. 












48 Hill a=2 



C3 89 Barlat a=8 
0 79 Hosford a=8 



Material 

FIG. 8. PREDICTED DEPENDENCE OF EARING (E45 
AND E90) ON THE CHOSEN YIELD CRITERIA FOR 
DIFFERENT A PIP. 





TABLE 3. CALCULATED STRENGTH RATIOS FOR 




Using the same "flange in uniaxial compression’' 
intuition that works qualitatively for E45, we can understand or 
at least explain why E90>0 for all three criteria; even though 
Y/X> 1 giving less hoop compressive strain in the 0 degree ear, 
P= 3 so there is much more lateral expansion, in the radial or ear 
direction for the 0 degree ear, as compared with die 90 degree 
ear where we have an R=1 compression test. Further, we can 
see why B89 would show more E90, since it predicts a 
stronger 45 degree hoop strength W/X. Less hoop compression 
at 45 degrees means more at 0 and 90 degrees and an 
exaggerated E90 compared with H48. This logic cannot explain 
the much-reduced E90 predicted by H79 with DYNA3D. In 
fact, since H79 predicts a lower F/X== 1.052, we would expect 
more hoop compression at 0 degrees and thus a higher value of 
E90 and perhaps even B89. A simple earing program was 
written to show this, and the results are shown in Table 4: 



TABLE 4. TRENDS IN E45 AND E90 EARS PREDICTED 
BY UNIAXIAL COMPRESSION EAR ANALYSIS. 








E90= 



H48 


H79 


B89 


.414 


.238 


.295 


-.080 


-.034 


■BQgQII 


.040 


.117 


.0451 



In essence, this simple earing program envisions the 0, 
45, and 90 degree directions in a cup undergoing a uniaxial 
compression test. The 45 degree direction is weighted twice, 
and we then demand a hoop force balance as the rim of the 
blank as a whole must see a hoop strain of e H = In (d p / d b ), 
where <ip=punch diameter, and d6=blank diameter. The 
individual hoop strains £^, £^ 5 , and £^° are adjusted to 
achieve a force balance in the hoop direction; as the magnitude 
of each of the compressive e d H increases, strain hardening 
increases the respective a e H . Further, thickening (depending on 
the local R(0+9O) value, as 1/(F(0+9O)+1)), and widening (as 
R(0-f9O)/(R(0+9O)+l)) increases the hoop force F®. The 
width and thickening effects can be modulated by adjustables 
a w and a t , where 0 < cq < 1 . As shown in Table 4, with 
a^—0.5 and a, = 1.0 as a best fit, this very simple ear 
analysis is able to reproduce, with fair quantitative agreement, 
the trends in E45 as shown by DYNA3D. The simple analysis 
also predicts that B89 should indeed have higher E90 than 
H48. The "higher” B89 dependence now seems to make sense. 
However, regardless of the values of a w and a t chosen, the 
simple ear analysis says that H79 should show the highest E90 
of all; this is the opposite of the finite element predictions. And 
yet, examining the strength ratios in Table 3, it does seem that 
E90 could be greater for H79, since the weak ear at 0 degrees 
will be squeezed out taller for H79 versus H48 or perhaps even 
B89. Either DYNA3D is leading us astray or a more complex 
ear analysis is required to capture the E90 trends. Such an 
analysis has been presented in {Van Houtte, 1993); a simplified 
extension of the above is added here in that spirit. In simple 
terms, we assume that although part of the deforming cup is in 
uniaxial compression, much of it has at least some component 
of uniaxial tension in the radial direction. In fact, there is some 
component, picked by another adjustable <x R in this analysis, of 
uniaxial tension all the way out to the flange. Assuming that 
material is flowing in the respective 0, 45, and 90 degree 
directions, we can find the uniaxial tension flow stress in the 
radial direction as a function of W/X or Y/X from 
eqns. (1 1-13); this is further modified by strain hardening due 
to plastic work assumed to occur in the hoop direction at these 
respective locations. (Naturally a coupled analysis would be 
more accurate, and perhaps require less adjustables, but 
separating the terms and effects does allow an increased 

intuitive feel). We can then back out a pseudo-radial strain 

from this pseudo-radial stress a R . The result of this, as shown 
in Table 5, is a strong effect on E45 and E90; die correction to 
E90 is such that for any value of adjustable a R , ( a R =0.25 
used in Table 5) the trends of E90 now agree with DYNA3D; 
that is, H79 shows the least E90 of all. 



TABLE 5. TRENDS IN E45 AND E90 EARS PREDICTED 
BY UNIAXIAL COMPRESSION PLUS RADIAL TENSION 
EAR ANALYSIS. 




Clearly, our simple intuition may fail us as the nature of 
earing and yield criteria become more complex. However, it 
seems worthwhile to find at least a qualitative simple solution 



















that allows our intuition to explain the finite-element or 
experimental trends. 



CONCLUSIONS 

Numerical simulations of the Limiting Draw Ratio and 
Earing phenomena during cupping have shown that many of the 
trends for fully "cross-rolled" materials, i.e. R = P , are 
intuitive, whereas the trends created when R & P (more typical 
or real sheet materials) are often not intuitive. In general, the 
H79 criterion seems to correlate best with available earing data; 
B89 is also reasonable and has a relatively simpler 
implementation. Both of these criterion give good correlation 
with data on LDR dependence. More work is needed on 
experimental and numerical trends of earing when R * P to see 
what is the best fit to trends of the data, in light of the complex 
earing interactions that occur in this circumstance. 
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